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Computing the viscosity of the QGP on the lattice 

Harvey B. Meyer 1 '*) 

1 Center for Theoretical Physics, Massachusetts Institute of Technology 
. 77 Massachusetts Avenue, 02139 Cambridge, USA 

o : 

f*"^ ■ I review the recent progress made in calculating shear and bulk viscosity on the lattice, 

£SJ ' and discuss ways to improve the calculation. 

ki2 ■ §1. Introduction 

S: 

' Models treating the system produced in heavy ion collisions at RHIC as an ideal 

fluid have had significant success in describing the observed large elliptic flow. 1 '' 2 ' 
i—j! Subsequently the leading corrections due to a finite shear viscosity r\ were estimated 3 ' 

and led to a remarkably small upper bound on -q. Recent relativistic viscous hydrody- 
r " ^ ■ namics simulations 4 ' -6 ' have started to achieve a tighter control over the systematic 

sources of uncertainty on the extraction of the shear viscosity to entropy density 
ratio, r]/s. On the theory side, it is therefore important to compute the QCD shear 
viscosity from first principles to complete the picture. Furthermore, since the heavy 
\ ion collision program at LHC will probe the quark-gluon plasma at temperatures 

£> ' about a factor two higher, it is crucial to predict the shear viscosity at ~ 3T C , and 

from there to predict the size of elliptic flow, before experimental data is available. 
| A small transport coefficient is a signature of strong interactions; strong inter- 

actions in turn require non-perturbative computational techniques. In this talk I 
. present a lattice calculation of the thermal correlators of the energy-momentum ten- 

sor in the Euclidean SU(3) pure gauge theory, and discuss methods to extract the 
shear and bulk viscosity from them. Computationally, the calculation is challenging 
enough without the inclusion of dynamical quarks, and physically, the thermody- 
namic properties of the QGP not too close to T c do not depend sensitively on the 
flavor content. 7 ' In perturbation theory, 8 ' the ratio of shear viscosity to entropy 
density of the pure gauge theory only differs by about 30% from full QCD 9 ' at a 
fixed value of a s (rj/s is smaller in the pure gauge theory). It should be appreciated 
that this is not much for a quantity which is infinite in the Stefan-Boltzmann limit, 
and the difference is actually reduced when comparing the pure gauge theory and 
full QCD at a common value of T/T c . 

Lattice calculations of the shear viscosity 10 ' -12 ' are based on Kubo formulas 
(see e.g. 13 ' and 14 '), which relate each transport coefficient to the small frequency 
behavior of the retarded correlator. By analytic continuation (see Appendix A), 
the latter is related to the Euclidean two-point function of the conserved current. 
The significant progress made in 12 ' was to obtain the energy-momentum correlators 
with an accuracy of about 3% on (isotropic) lattices with temporal resolution up to 
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N T = 12. This progress was mainly due to the use of a multi-level algorithm. 15 ) 

Besides technical improvements that concern the discretization of the correlators, 
I present new ideas to enhance the sensitivity of the lattice data to the low-frequency 
region of the spectral function, which determines the transport properties. Subtract- 
ing correlators with different spatial momentum p or different temperature T will 
cancel off most of the contribution of the high-frequency modes. The task is then 
to solve for the difference of the spectral functions for two different values of p or 
T. This difference is no longer positive definite and solving for the spectral function 
requires new methods, such as developed in Ref. 16 ) 

The correlators computed on the lattice are (L$ = 1/T) 

C s (x ,p) = Llf d 3 x e*P' x (T 12 (0)T 12 (x ,x)}, 

C b (x ,p) = ^-Yl ^ (Tu(0)T n (x , X )). (1-1) 

i,j=i J 

For our purposes the spectral functions are then defined by 

coshw(iLo — xq 



f°° coshu;(iLo — xq) 
('a,b(x ,p) = L 0j Ps,b(u,p) s - mh ^Lo ( 1-2 ) 



The shear and bulk viscosities are given by 10 ) 

The spectral functions are positive, p(u,p)/u > 0, and odd, p(— u, 0) = —p(u,0). 
If not specified, p is set to zero in this talk. In Ref. 16), I defined the following 
moments of the spectral function (n = 0, 1, . . . ): 



(co 2n ) = I*[ 
Jo 



oo , ,2n n (, ,\ 

(1-4) 



. w*p(w) d 2n C 
diu — 



o sinhu;Lo/2 cIxq 



2n 



xq=L /2 



The latter equality implies that they are directly accessible to lattice calculations. 
A first observation is that in infinite spatial volume and in the continuum limit, 

c s ( Xo ,p) = ^f J d 3 xe*p x ((rn-r 22 )(o)(rii-r 22 )(xo,x)), P = (o,o,p). (i-5) 

All the results concerning C s presented in this talk have been obtained by discretizing 
this form. Secondly, because Tu = — Too, an d because J <i 3 x(Too(x)0(O)) c = 
T 2 8t{0)t for any local operator O and xq / 0, 

J d 3 x (T^(x) T vv (0)) c = T 2 8 T (e - 6P) + J ci 3 x (T u (x) T kk (0)) c . (1-6) 

We have used (Th)t-o = — 3P, (Too)r_o = £■ It is therefore convenient to study the 
two-point function Cq of the trace anomaly 9 = T^. 
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§2. Thermal correlators from isotropic lattices 

Figures (1,2) show the p = scalar and tensor correlators for the range of tem- 
peratures T c to 3.2T C . This data was obtained on N T = 8 lattices with the Wilson 
action and the 'bare-plaquette' discretization 17 ) of the energy-momentum tensor, 
and I implemented the tree- level improvement. 12 ) The leading cutoff effects are thus 
O(g 2 a 2 )- 

The tensor correlator exhibits near-conformal behavior, while in the scalar case 
large departures from conformality are seen, particularly near T c . The curves are 
the leading-order perturbative results. For the scalar correlator, 16 ) which is 0(a 2 ), 
a choice has to be made for the value of the coupling. The value that matches the 
LO prediction for e — 3P with the non-perturbative value 18 ) at 3.22T C is a s (2irT) = 
a* = 0.289. 16 ) I then use the one-loop evolution on Fig. (2), ci s (2itxq 1 ) = a*/(l — 
Tj^otg log(Txo))- The non-perturbative correlator is somewhat flatter than the LO 
perturbative prediction. A study of cutoff effects 16 ) for T c < T < 2T C shows that the 
N T = 8 data at xq = Lq/2 is an overestimate, at most by 20%, of the continuum 
correlator. Large finite-size effects can be excluded on the basis of the large aspect 
ratio, LT = 6. 




Figure (3) displays the first two moments, Eq. (jl-4|) . of the tensor spectral func- 
tion at p = 0. They are normalized by the leading perturbative result. 12 ) I also show 
the corresponding ratio of moments in the (conformal) AT = 4 SYM theory, 13 )' 19 ) 
computed by AdS/CFT methods. It is remarkable that the typical size and the sign 
of the deviations from the free approximation is the same in both theories. 

Based on the data displayed on Fig. (1,2) and data at smaller spacing, I es- 
timated 12 )' 16 ) the shear and bulk viscosity by expanding the spectral function lin- 
early in a set of orthogonal functions U((u). In this way I obtained for instance 
rj/s = 0.13(3) at 1.65T C . Due to the small size iVp of the set, the functions fail to 
satisfy the completeness relation by an amount quantified by the resolution function 
6(uj,uj') = Yl^=i u zd jJ ) u £d jJ ')- Figure (4) shows the resolution function for u//T = 0, 
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10 and 20 for N p = 4. It is broad, and quite far from resembling a delta function 
at <jj = u/. Nevertheless, in cases where the spectral function is smooth, such as the 
strongly coupled SYM theory, 13 )' 19 ) the method works well, see Fig. (4). 

Figure (3) demonstrates that even at 1.6T C , the deviation of the first moment, 
to which the transport peak contributes, from the non-interacting approximation, is 
only about 10%. This observation had previously been made in the finite-temperature 
perturbation theory framework 20 - 1 and in the strongly coupled SYM theory. 13 ) It 
therefore appears necessary to investigate the analytic structure of the spectral func- 
tion to understand where this lack of sensitivity to the low-energy degrees of freedom 
comes from. 
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Fig. 3. The lowest two moments of the ten- 
sor spectral function for N T — 8 on the 
isotropic lattice. 



Mock reconstruction of the N=4 SYM spectral function 




co/T 



Fig. 4. Reconstruction of the SYM spectral 
function from the N T = 8 Euclidean cor- 
relator; resolution functions at 3 points are 
shown. 



§3. Analytic structure of the spectral function 

The leading-order perturbative behavior of the spectral functions can be found in 
Ref. 12), 16), they grow as w 4 at high frequencies, and at T = the NLO corrections 
in a s are known. 21 )' 22 ) Furthermore, at T = 0, Lorentz invariance implies 

p(u, p) = sgn(cu)0(cu 2 - p 2 ) p{y/u* - P 2 , 0) (3-1) 

in the scalar channel. Based on the operator product expansion, 23 ) we expect the 
perturbative series of the w 4 coefficient to be independent of temperature. Therefore, 
this term, which makes a large contribution to the Euclidean correlators without 
telling us anything about thermal physics, can be eliminated by subtracting the 
T = spectral function from the finite-temperature one. 

For low momenta and frequencies, hydrodynamics predicts the functional form 
of the spectral functions in the shear channel and the sound channel {pu n, defined 
as p s with T\2 replaced by Txi). For p = (p, 0, 0), and v s being the velocity of sound, 

Ps{u,p) V 

u itlo 2 + ( V p 2 /(Ts)) 2 ' 1 ' 
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(3-3) 



Vs V) 2 + (wp 2 (|?? + C)/(rs)) 2 ' 



See Ref. 13) for a particularly clear derivation. It is therefore of interest to study also 
correlators with non-vanishing spatial momentum. Ultimately, observing this struc- 
ture in the spectral function is the best way to give us confidence in the extraction 
of the viscosities. 

Finally, we remark on a subtlety in the calculation of bulk viscosity. On general 
grounds, pb is expected to not have any delta function at to = in an interacting 
theory. This would indeed reflect the conservation of (part of) the momentum cur- 
rent, which would imply in particular that such a current never dissipates and the 
system never reaches equilibrium. The spectral function for the Cg correlator, in 
view of Eq. (|l-6p . must then contain the term T8t(£ — 6P)uj5(uj). This singular term 
was missed in. 16 ) 

§4. Methods to enhance the sensitivity to the low-frequency region. 

4.1. Strategy I: exploiting the T = spectral function 

The idea is to solve the integral equation fjl • 2|) for the linear combination 



In words, subtract the zero-temperature spectral function, except for its one-particle 
contributions. Indeed there are two glueballs below the two-particle threshold in 
both the scalar and the tensor channel. In infinite spatial volume, the function 
subtracted from p(T,u>) is exactly zero below 2M ++ ~ 3GeV. The low-frequency 
region is therefore unaffected, but the high-frequency asymptotics of the function to 
be reconstructed is now u>°, up to logarithms. This is a dramatic improvement. 

The first step is thus to determine the T = spectral function, on which Lorentz 
symmetry places much stronger constraints. The Euclidean correlators deep in the 
confined phase are shown on Fig. (5). They are computed with the isotropic Wilson 
gauge action, with a ~ 0.068 and 0.051 fm. The correlators fall off rapidly at large 
distance, where a significant signal is obtained only due to the multi-level algorithm. 
The correlators have been treelevel improved. 12 )' 16 ) After this improvement, the two- 
loop perturbative prediction 21 ) of the tensor correlator and the one-loop prediction 
for the scalar channel (both with a s set to 0.25) compare rather well with the data 
at short distance. At large distance xq > 0.5 fm, the data is compared to the contri- 
bution from pip, namely that of the two stable glueballs present in each channel. I 
took their masses from 25 ) and computed their matrix elements separately. 27 ) 

4.2. Strategy II: linear combinations o/p/0 spectral functions 

From Eq. (|3-ip . we expect that for u 3> T, the u> 4 contribution cancels up to 
two-loop order in the linear combination 



Ap(T,uj) = p(T,u) - 0(0,a;) - pi v (u)). 



(4-1) 



-e=) + ^^p), 6>l. 



(4-2) 
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In the tensor channel, this means that the leading large-u; behaviour of this linear 
combination is 0(a 2 w 4 ). Figure (6) displays linear combination f|4-2j) of the cor- 
relators C s for 6 = 4 and p = (0, 0, irT). A cancellation by almost two orders of 
magnitude takes place and the data is consistent with zero at all xq. This is a re- 
markable fact; by contrast C s (xn,0) — C s (xo, ^p) does not vanish. Obviously such 
linear combinations deserve further investigation. The data on Fig. (6) was obtained 
on an anisotropic lattice with the clover discretization 17 ) of the energy-momnentum 
tensor, with non-perturbative normalization factors determined using the thermo- 
dynamics data in Ref. 26). 
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Treelevel Improved Correlators at T=T c /2 (20 4 , 28 4 ) 
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Fig. 5. C s and Ce in the confined phase. 
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Tensor correlators (10x20 , %=2, clover discretization) 

C s (x n , 0) - (4/3) C s (x , p) +(l/3)C s (x , 2p) 

C s (x , 0) - C s (x () , p) — ^- 
C s (x , 0) / 100 — • — 



0.2 



T=1.49T C LT=4 
p=(0,0,p 3 ), p 3 = kV2 



0.3 0.4 

Tx„ 



0.5 



Fig. 6. Treelevel improved tensor correlators. 



§5. Conclusion 

The correlators of the energy-momentum tensor contain information on the 
plasma that is complementary to that of thermodynamics. They are sensitive to 
the low-energy degrees of freedom rather than the bulk of them, in particular they 
contain the information on the transport properties of the system. 

Describing the Euclidean data with a smooth spectral function consistent with 
positivity, parity and the perturbative large-frequency prediction leads to a low shear 
viscosity to entropy ratio, 12 ) 1 < Anri/s < 2 in the temperature range 1.2 < T/T c < 
1.7. 

In order to increase the sensitivity of the Euclidean correlators to the low- 
frequency domain of the spectral function, I have proposed two different strategies 
that aim at subtracting the contributions of high-frequency modes. The goal is to 
challenge the smoothness assumption made on the spectral function. 12 )' 16 - 1 

A central question is whether perturbation theory can explain most properties 
of the plasma at ~ 3T C , a typical temperature probed at LHC. If it is to account 
correctly for the viscosities, then it should accurately describe the corresponding 
Euclidean correlators. In the scalar channel at 3.2T C , the agreement of LO perturba- 
tion theory with the lattice data depends sensitively on the choice of coupling value. 
Figure (3) suggest that the agreement is good in the tensor channel at 3.2T C , how- 
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ever the data is also compatible with a strongly coupled scenario, as the qualitative 
comparison with AdS/CFT reveals. The methods I presented to subtract the UV- 
contributions have the potential to elucidate which is the more appropriate picture 
of the plasma at temperatures typical of the LHC experiments. 
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Appendix A 

Analytic continuation of the retarded correlator 

Our goal in this section is to relate the spectral function of a conserved opera- 
tor, defined via a Euclidean correlator, to the imaginary part of the retarded Green's 
function in frequency space. A Kubo formula relates the latter to a transport coef- 
ficient of the finite-temperature system, in the case of a conserved operator; for the 
shear viscosity, r\ = — lim^-^o ^ Im G 1 r 2,12 (uj), where G R 2 ' 12 is the retarded Green's 
function of T 12 (see Ref. 13), 14)). 

The Euclidean correlator, Cs(t) = (O(t)O(0)}, t > 0, has the spectral represen- 
tation 

Ce(x ) = I E \Onm\ 2 e- LoE -e E - x °. (A-l) 

n,m 

Here O nm = {n\0\m), E nm = E n — E m . One easily finds that C^(t) can be expressed 
in terms of the spectral function p(Lq,uj) (see Eq. II -2 j) 

P(L ,lj) = ^smh(u;Lo/2)^5(u - E nm )e-^ +E ^ 2 \O nm \ 2 . (A-2) 

n,m 

On the Minkovsky side, the retarded correlator iG R (t) = 6(t) ([0(t), O(0)]) has 
the spectral representation 

iG R (t) = ^ \O nm \ 2 e- LoE " {e iEnmt - e ~ iEnmt ). (A-3) 

n,m 

It is related to the Euclidean correlator by the relation 

iG R {t) = lim{C E (it + e) -C E (-it + e)), t > 0. (A-4) 

A small positive real part in the argument of Ce guarantees the finiteness of the 
expression. In terms of the spectral function, we obtain 

/oo 
du p(L ,u)e- iuJt e-^ e , t > 0. (A-5) 
-oo 
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We have exploited the fact that p is odd in uj, a property manifest in Eq. (|A-2I) . 
The Fourier transform of Gr, Gr(lo) = f^°dt e lujt Gn(t), converges if we give its 
argument a positive imaginary part: 

Gr(o; + zo) = - / do; — , a; real. (A-6) 

J-oo u? -u-%6 

In particular, 

/OO / 1 5 
da;' (-7r/)(L ,w') e" |w ' |e ) - — N2 ■ x2 = -np(L ,u), 
vr (a;' - a;)^ + ^ 

(A-7) 

where we have recognized one of the standard representations of the delta function, 
and let e — > in the last step. 
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